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Diffusion-limited aggregation has a natural generalization to the "77- models" , in which 77 random 
walkers must arrive at a point on the cluster surface in order for growth to occur. It has recently 
been proposed that in spatial dimensionality d — 2, there is an upper critical 7/ c = 4 above which 
the fractal dimensionality of the clusters is D = 1. I compute the first order correction to D for 
r] < 4, obtaining D = 1 + |(4 — rj). The methods used can also determine multifractal dimensions 
to first order in 4 — r\. 
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The formation of patterns in nature is often controlled by diffusive phenomena. The branching of a less viscous fluid 
such as water injected into a more viscous fluid such as oil, the dendritic complexity of a snowflake, or the formation 
of veins of metals on the surface of certain rocks, all display pattern formation processes controlled by the diffusive 
transport of some quantity |jj . 

The simplest model for diffusion-controlled growth was introduced twenty years ago by Witten and Sander Q . Their 
model of "diffusion-limited aggregation" (DLA) describes the formation of an aggregate by sequential deposition of 
randomly walking particles arriving from infinity. There is an electrostatic formulation of this process, in which the 
n-particle cluster is chosen as an equipotential of a Laplacian field, which has a source at infinity. The local growth 
probability, i.e., the probability of deposition of the n + l'st particle, is then chosen proportional to the local electric 
field on the surface of the cluster. 

A natural generalization of this model is to fix the growth probability on the surface proportional to the 77'th power 
of the local electric field. This corresponds to a random walk model in which 77 independent walkers must arrive at 
the same surface point in order for growth to occur at that point. These "77-models" were originally introduced by 
Niemeyer et al. as models for dielectric breakdown, and represent a useful formal extension of DLA ||. 

These models were used in an important recent work of Hastings to propose a systematic perturbative approach 
to DLA |Q. This work argued that the fractal dimension D of 77- model clusters collapses to D = 1 in spatial 
dimensionality d = 2 for 77 > 4, and that this value of rj c = 4 therefore represents an upper critical 77 for these models. 
Dimensions and other properties of models for rj < 4 can then be determined by perturbative renormalization in 
4 — rj. The case of DLA (77 = 1) is in principle accessible, although satisfactory agreement with the numerical result 
D = 1.71 may be difficult to achieve given the large value of 4 — 77 required. However, considerable computational 
difficulties arose in implementing this program. Nevertheless, rough numerical results for the first order correction to 
D were obtained, which agree with the exact result expressed in Eq. (|l|) below. 

Many of the ideas used by Hastings originated in the "branched growth model," a phenomenological treatment of 
DLA proposed a number of years ago by my collaborators and myself ]s| . The purpose of the current work is to show 
that the branched growth model actually allows easy computation of perturbative terms, at least to first order in 
4 — rj. This ease can be understood as a consequence of the branched growth model becoming exact as one approaches 
rj = 4. In particular, 1 obtain the result that the dimension D is given to first order in 4 — 77 > by 

# = l + i(4 - 77)+ 0(4 - 77) 2 . (1) 

1 obtain as well first order expressions for the multifractal dimensions of the growth measure. 

In this work, I will first review the salient features of the branched growth model, and I will argue that to lowest 
order in 4 — 77 it well represents the dynamics of the underlying Laplacian growth process. I compute the dimensions 
of the clusters by two different means, and show that both methods give Eq. (Q). I then derive an integral formula 
for the multifractal dimensions to 0(4 — rj), and give a simple approximation to the actual values of these dimensions 
at this order. Finally, I discuss prospects for extending this computation to higher orders in 4 — 77. 

The branched growth model places a fundamental importance on the microscopic process of tip-splitting, whereby a 
growing branch forks into two growing branches. This process occurs at a microscopic scale, on the order of the particle 
size or cutoff a in dimension. Thus the frequency and detailed dynamics of tip-splitting is controlled by microscopic 
and presumably non-universal details of the way particles attach at or near the tip of a growing branch. We regard 
tip-splitting as the fundamental stochastic process in the model; we disregard all other forms of stochasticity such as 
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the "shot noise," i.e., the purely statistical variations in the number of particles depositing at different positions in 
the cluster. The reader should note that the precise role of stochasticity in DLA has recently been quite controversial 
||; although we believe that the theory to be outlined in this work is stable against obvious additional sources 
of stochasticity such as shot noise, a complete understanding of the roles of different kinds of noise requires more 
systematic study. 

Once a branch splits into two, we follow its additional development by implementing a deterministic version of the 
?7-model growth rule. Near the tip of a linear equipotential, the electric field in two dimensions diverges as 

E(w) ~ w~ 1/2 (2) 

with w the distance from the tip. Thus, for 77 > 2 the growth measure, which is proportional to E rt {w), is entirely 
dominated by growth at the tips. Since we work near 77 = 4, we need only follow the progress of the tips in the 
deterministic portions of the growth, as well as keeping track of the generation of new tips through (stochastic) 
tip-splitting. 

Consider two branches emanating from the same tip-splitting event. The masses (particle numbers) of the two 
branches we will write as til for the left-hand branch, and tir for the right-hand branch; the growth measures of the 
left-hand and right-hand branches we write as pl, Pr respectively. Defining relative growth rate and mass parameters 
x and y respectively by 

(3) 
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with tit — til + tir, we see by elementary manipulations that 

dy 



dlog(riT 

In principle, we can define a function G of the overall cluster geometry such that 

dx 



y- (5) 



d\og(n T ) 



= G. (6) 



In the branched growth model, G is taken to be a function of x and y alone. 

Let us consider a growing fork, i.e., a branch with equal sub-branches (the tines of the fork) immediately after 
tip-splitting. It is convenient to describe the fork by the conformal map w = F(z) that maps the real axis in the 
z-plane onto the fork in the physical w-plane; \dwfdz\~ 1 gives the local electric field at w. If we choose a fork for 
which the angle between the tines is 9%, the map is given by 

w = F(z) = z ai (z 2 - 1) Q2 (7) 

with ai = 9i/tt and ai = 1 — (61 /2w), where the constant 1 is chosen arbitrarily. The derivative of the map F'{z) 
possesses zeroes at z± — ±a/oi/2, which correspond to the points of the fork. We can fix cei by requiring that the 
points are oriented towards the maximum field, so that the fork geometry is unchanged by the growth process; this 
requires that 



d ( F'(z) \ d ( F'{z) 



dz \z — z+) J dz \z — z-) 



(8) 



which after some algebra determines ai = | or 6% = g7T 

The competition between the two growing tines of the fork is intrinsically unstable. A simple computation shows 
that the eigenvalue of the instability is j| 

" = f - I- (9) 

Since the two tines are supposed to be created with approximately constant probability to be found near the fixed 
point (x,y) = (i, 5), the eigenvalue of the instability can be related to the probability that the branch pair is still 
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active (i.e., one branch has not been entirely screened by the other) as t%t grows; this probability is P oc rii^ v . Since for 
a main branch of length r there will have been 0(r) opportunities to branch, requiring that 0(1) sidebranch always 
be active implies v = D^ 1 Q, which already suggests Eq. (|l|). However, it is productive to consider this question in 
more detail. 




FIG. 1. The geometry of the growing branch. After each tip-splitting event, the weaker branch is screened by the stronger. 
The tip-splittings are indexed by j > 0. Tip-splittings are separated by m particles on the main branch. 

Consider a growing branch of n particles, which tip-splits every m <C n particles. The side-branches thus generated 
persist for a certain distance, and are then screened by the main branch. Index these tip-splittings by < j < J 
(see Figure 1). Then at each j there can be defined parameters Xj, yj, and rij_T, giving the relative mass and growth 
probabilities of the side-branches, as well as the total mass of the remainder of the main branch plus the side-branch 
in question (the total mass to the right of the branch point in Figure 1). We choose our definition of "left" and "right" 
in Eqs. (^j^) so that each Xj < |. The growth rate at the overall tip is then given by 

j 

Pti P = l[(l-x j ). (10) 

3=0 

Suppose that at the j'th branching, the initial value of xj is given by 

^i,T = l) = ^-^ (11) 

defining a random variable for the branching Ej. The distribution of e, p(e), is chosen so that Xj does not have any 
singularity in its initial distribution near x = 1/2-this reflects the microscopic origin of the stochasticity. 

]imp{s)=p s v -\ (12) 

with the large e behavior constrained by the integrability of p(e). In the branched growth model, the dependence of 
Xj on rij^T and ej at each branching point is the same; we now assume this to be the case in this more general model 
as well. The argument below will establish that this assumption is correct to lowest order in 4 — 77. 

If the dynamics of each branching point are the same, then there must exist some G(x,y) in Eq. (^]). The values 
of x and y can then be integrated to obtain x(h), y(h) from Eqs. (|^-^|), with h oc tit and (x(0),y(0)) = (|, \). Note 
that this choice of variables means that the dependence of x on e can be encoded by 

x(riT,e) ~ x(en,T = h), (13) 

and similarly for y. This formula becomes exact for nx large. 
We can now see that to O(po) 

(pup) = 1 - 53<aBj> = 1 - / dzpoe v - l x{en jtT ), (14) 
3 3=0 J 

or, extracting the dependence on {tIj.t} from the integrals, 
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EKt> {p°j™ dh h "~ Xx ( h ) 



(Pup) = 1 - 5Z(«7t) A) / * h"- x x{h) . (15) 



To evaluate X^=o( n j t)' we use a smT pl e trick. Let us first assume that 

EKr) = Alogn. (16) 

We will now confirm this form, and compute the parameter A. Consider the second branching at j = 1. We have 

since the number of particles in the main branch below the second sidebranching is n\ — (1 — yo)n — m (see Figure 
1). Note that n = Uqt- This leads to 

— = -^A / d/i ^- 1 log(l - j/W) + — <(1 - yo)- 1 )- (18) 

As /i — > oo, the weaker branch will die at some fixed mass and y(h) — > h/h. Thus as v — > 1, the integral over /i on 
the right-hand side of Eq. ([l8]) diverges as 

dhh"- 1 log(l-y(h)) = -j^+0(l) (19) 
Thus we obtain immediately 

A = IzJL + 0(1 - v f (20) 

so that the sum in Eq. (Esh is 0(1 - z/) = 0(4 - /?). 

The fact that this "propagator" sum is oc (1 — ^)logn is the key formal result of this work, which allows us to 
construct a direct renormalization group for the dimension and other properties of the 77-models. The procedure, in 
principle, is as follows. First, a naive perturbative expansion in po is constructed, along the lines of ref. The 
computation above shows an example to first order in pq. This expansion should account both for the different 
contributions of the various tips to the quantity being computed, as well as the influence of the internal structure 
of the various branches on the functions x(h), y(h). In this expansion, sums over {nj!^) such as that appearing in 
Eq. (^) will appear, as well as more complex, albeit still ultimately logarithmic, sums. Performing these sums, one 
will replace the original series in pq with a logarithmic series in 1 — v. The methods of ref. § easily show that the 
higher order terms in pg will be higher order in 1 — v upon computation of these sums. This series then forms the 
basis of a direct renormalization calculation of the quantity of interest. 

We illustrate this by returning to the growth rate of the overall branch tip. To zcroth order in 1 — is, we can represent 
the structure of the two branches, which we use to compute x(h) and y(h), by a fork with two growth sites at the 
fork tips. Thus we generalize the fork with equal-length tines to a fork which may grow different lengths on the two 
sides-in so doing, there will be a tendency for the two sides to curve away from being perfectly straight and separated 
by an angle 9\, as shown in Figure 1. Standard techniques of integrating the conformal map for Laplacian growth 
structures will suffice for determining the full x{h) and y(h) for this case M; below we give a simple approximation 
to these quantities. To this order, there are no additional variables describing the internal structure of the branches, 
so we are justified in assuming the same xQi) and y(h) at each branch point. 

However, it turns out that we need not compute x(h) or y(h) explicitly in order to compute the first order correction 
to the tip growth rate. From the definitions of x(h) and y(h), we have that 

z, (21) 

driT 

implying by integration of parts that 

o />oo 

dh h v ~ l x{h) = (1 - v) I dh h u ~ l y{h). (22) 

n JO 
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The divergence of f dhh lJ 1 log(l — y(h)) as v — > 1 originates from large values of h, or small values of y. Thus to 
lowest order in 1 — v we have 



oo 



and 



so that 



dh h u - l \og{l - y(h)) = - / dhh^yQi), (23) 



f°° On f°° dh h 1/ ~ 1 x(h) 

5>7>o / dft h^x(h) = ,00 ,„ i , logn = (1 - i/Jlogn, (24) 
Y ^° ^°io dhh v ^(h) 



p Up = l-(l-v)\ogn + 0(l-vf. (25) 
Since we expect p^ p oc n D with Z) mass-radius scaling dimension Jicj ], this implies with Eq. (^) that 

D = 1 + 1(4-7,) + 0(4 -t?) 2 , (26) 

as advertised. 

To compute multifractal dimensions to 0(4 — rj), I use similar techniques. Following ref. Q, and using z as an index 
to growth tips, we see that the multifractal spectrum cr(q) for the growth measure (not the harmonic measure) is 
given by 

(%2p1) sn"^ = l + Alognj^ ctt + (1 - x(h))« - 1]| + 0(1 - ^) 2 , (27) 



a(g) = ^ / ^[^(/i) + (l-x(/i))«-l]^+0(l-^) 2 . (28) 



yielding 



To obtain explicit results for the multifractal dimensions, we need the trajectories x(h) and y(h). I use a simple 
artifice that gives a useful approximation to this trajectory. The integral in Eq. ( |28] ) is dominated by values of x 
near x = \. We can thus approximate the integral by taking the linear trajectory in the x — y plane near the center 
(x, y) = (|-, h) and extending it to the boundaries x = 0, 1. (This approximation is referred to as "model Z" in ref. 0.) 
Explicitly, we write, taking the lowest order v = 1, 

f - — h ifh< - 

^Ho, S£>|, (29) 



and 



^={i/ 8 f' !f htt ( 3 °) 

which agrees with Eq. (|). 

The approximate result for the multifractal dimensions is then 

o(q) = 2(4-7?) ^Y + °(4-^) 2 - ( 31 ) 

Note that due to our total suppression of non-growing portions of the measure, we do not recover the identity 
cr(0) = —D, 

The extension of these results to higher orders in 4 — 77, and in particular to the case of DLA (77 = 1), will require 
some further formal development. Reference |j) successfully computes the most divergent and next-most divergent 
terms at all orders in po for the multifractal dimensions for the branched growth model; the behavior of the higher 
order logarithms in this case does allow resummation of the theory for, e.g., quenched and annealed multifractal 
dimensions. In our case, we need to add a family of terms representing the deviations from the perfect branched 
growth model behavior, which arise from fluctuations in the internal structure of the branches. Fortunately, there are 
indications that these fluctuations are also renormalizable || . 

I am grateful to M.B. Hastings for drawing my attention to ref. 0], and for a critical reading of this manuscript. 
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